The End of Unified Dark Matter? 
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Despite the interest in dark matter and dark energy, it has never been shown that they are in fact 
two separate substances. We provide the first strong evidence that they are separate by ruling out a 
broad class of so-called unified dark matter models that have attracted much recent interest. We find 
that they produce oscillations or exponential blowup of the matter power spectrum inconsistent with 
observation. For the particular case of generalized Chaplygin gas models, 99.999% of the previously 
allowed parameter space is excluded, leaving essentially only the standard ACDM limit allowed. 
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I. INTRODUCTION 

Despite the broad interest in dark matter and dark en- 
ergy, their physical properties are still poorly understood. 
Indeed, it has never even been shown that the two are in 
fact two separate substances. The goal of this paper is 
to show that they are. 

There is strong evidence from a multitude of observa- 
tions that there is about six times more cold dark mat- 
ter (CDM) than baryons in the cosmic matter budget, 
making up of order 30% of critical density [1,2]. In ad- 
dition to this clustering dark component, observations of 
supernovae, the cosmic microwave background fluctua- 
tions and galaxy clustering provide mounting evidence of 
a uniformly distributed dark energy with negative pres- 
sure which has come to dominate the universe recently 
(at redshifts z < 1) and caused its expansion to accel- 
erate. It currently constitutes about two thirds of the 
critical density [1-3]. 

Although the dark energy can be explained by intro- 
ducing the cosmological constant (A) into general rela- 
tivity (lending the standard model the name ACDM), 
this "solution" has two severe problems, frequently trig- 
gering anthropic explanations and general unhappiness. 
The first problem is explaining its magnitude, since the- 
oretical predictions for A lie many orders of magnitude 
above the observed value. The second problem is the so- 
called cosmic coincidence problem: explaining why the 
three components of the universe (matter, radiation and 
A) presently are of similar magnitudes although they all 
scale differently with the Universe's expansion. 

As a response to these problems, much interest has 
been devoted to models with dynamical vacuum energy, 
so-called quintessence [4] . These models typically involve 
scalar fields with a particular class of potentials, allowing 
the vacuum energy to become dominant only recently. 
Although quintessence is the most studied candidate for 
the dark energy, it generally does not avoid fine tuning 
in explaining the cosmic coincidence problem. Recently 
several alternative models have been proposed, such as 



the condensate models of [5]. 

An alternative to quintessence which has attracted 
great interest lately is the so-called generalized Chap- 
lygin gas (hereafter GCG) [6,7,9,8,10,11] (see also the re- 
lated earlier work of [12]). Rather than fine tuning some 
potential, the model explains the acceleration of the Uni- 
verse via an exotic equation of state causing it to act like 
dark matter at high density and like dark energy at low 
density. The model is interesting for phenomenological 
reasons but can be motivated by a brane-world interpre- 
tation [8,7]. An attractive feature of the model is that it 
can explain both dark energy and dark matter in terms 
of a single component, and has therefore been referred 
to as unified dark matter (UDM) or "quartessence" [13]. 
(See also [17].) 

This approach has been thoroughly investigated for 
its impact on the 0th order cosmology, i.e., the cos- 
mic expansion history (quantified by the Hubble parame- 
ter -ff [z]) and corresponding spacetime-geometric observ- 
ables. An interesting range of models was found to be 
consistent with SN la data [13] and CMB peak locations 
[15]. 

Some work has also studied constraints from 1st order 
cosmology (the growth of linear perturbations), finding 
an interesting range of models to be consistent with cos- 
mic microwave background (CMB) measurements [16]. 
There is, however, a fatal flaw in UDM models that man- 
ifests itself only at recent times and on smaller (Galactic) 
scales and has therefore not been revealed by these stud- 
ies. As we will see, this flaw rules out all GCG models 
except those that are for all practical purposes identical 
to the usual ACDM model. 

The rest of this Letter is organized as follows. In the 
next section, we review the fundamentals of the GCG 
model. We then consider in section III the evolution of 
density inhomogeneities in the model and use the pre- 
dicted matter power spectrum to constrain it with ob- 
servational data. Finally, we describe how the basic flaw 
that rules out the GCG model is indeed a generic feature 
of a broad class of unified dark matter models. 



1 



II. THE CHAPLYGIN GAS 

A standard assumption in cosmology is that the pres- 
sure of a single substance is, at least in linear perturba- 
tion theory, uniquely determined by its density. A gener- 
alized Chaplygin gas [6,8,9] is simply a substance where 
this relation p{p) is a power law 
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with A a positive constant. The original Chaplygin gas 
had a = 1. The standard ACDM model has two separate 
dark components, both with a = — 1, giving a constant 
equation of state w = p/p that equals for dark matter 
and —1 for dark energy. 

By inserting equation (1) into the energy conservation 
law, one finds that the GCG density evolves as [13] 
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where a(t) is the cosmic scale factor normalized to unity 
today, i.e., a = (1 + z)^ 1 where z denotes redshift. Here 
B is an integration constant. The striking feature here 
is that although the GCG has p oc a~ 3 when sufficiently 
compressed, it's density will never drop below the value 
^i/i+a nQ ma ^ er now mu ch you expand it. Defining 
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p* = (A + B) — 
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equation (2) takes the form 

P (a) = p, [(l - n*j + n* m a-* 1+a A 



(3) 
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For comparison, a standard flat model with current CDM 
density parameter Q m as well as dark energy density (f — 
O m ) whose equation of state w* is constant gives 



p(a) = p* \{l - n m )a-^ 1+w -^ + fl m a- 3 



(5) 



We see that the last two equations bear a striking simi- 
larity even though the former involves a single substance 
and the latter involves two. Both have two free param- 
eters. Both have the current density p(l) = p*. Making 
the identification Sl^ = fl m , both have p(a) — > r2 m /9*a~ 3 
at early times as a — ► (for u>* < 0), showing that £1^ 
can be interpreted as an effective matter density in the 
GCG model. Indeed, for the special case a = and 
w* = —I, we see that both models coincide with standard 
ACDM. For a = the GCG model becomes equivalent 
to ACDM not only to 0th order in perturbation theory as 
above but to all orders, even in the nonlinear clustering 
regime. 

The 0th order cosmology determined by equation (4) 
together with the Friedman equation 
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(which determines a(t) and the spacetime metric to 0th 
order) has been thoroughly in previous work [6,8], and by 
studying constraints from supernovae observations, Mak- 
ler et al. [13] have placed interesting constraints on the 
(a, fijjj parameter space. 



III. GROWTH OF INHOMOGENEITIES 

Let us now consider the evolution of density perturba- 
tions in this UDM model. Following the standard cal- 
culations of [17], we obtain for the relativistic analog of 
the Newtonian 1st order perturbation equation in Fourier 
space that a density fluctuation 5k with wave vector k 
evolves as 



S k + H6 k [2-3(2w-c 2 s )] 
-^H 2 5 k [l - 6c 2 s - 3w 2 + 8w] 



(7) 
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where the equation of state w = p/p and the squared 
sound speed c 2 = dp/ dp are evaluated to 0th order and 
hence depend only on time, not on position. (We use 
units where the speed of light c = 1 throughout.) This 
equation is valid on subhorizon scales |k| ^> H/c. In 
other words, the growth of density fluctuations is com- 
pletely determined by the two functions w(a) and c 2 (a). 
Combining equation (1) and equation (4), these two func- 
tions are [13] 
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This shows a second reason why the GCG has been con- 
sidered promising for cosmology it starts out behaving 
like pressureless CDM (w k> 0, c s w 0) early on (for 
a«l) and gradually approaches cosmological constant 
behavior (w w —1) at late times. There is also an in- 
termediate state where the effective equation of state is 
p = ap [6]. (Going beyond 1st order perturbation theory, 
the GCG that gets gravitationally bound in galactic ha- 
los maintains its density high enough to keep acting like 
CDM forever.) 

To solve equation (7) numerically, we change the inde- 
pendent variable from t to In a. Using the properties 



dt d In a ' 
where ' = d/dlna and 
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equation (7) takes the form 
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the pressure term in equation (13) becomes simply 
(X c k) 2 , so we expect oscillations or exponential blowup 
in the power spectrum on scales k > A" 1 . These are 
created mainly during the recent transition period when 
both a and — w are or order unity (growing from to 1), 
and since neither effect is seen in observed data, we there- 
fore expect to obtain constraints of order \a\ < (H/k) 2 , 
the squared ratio of the perturbation scale to the horizon 
scale. This heuristic argument thus suggests that Galaxy 
clustering constraints on scales down to 10/i _1 Mpc would 
give the constraint \a\ < (lOh-^Mpc/SOOO/i^Mpc) 2 « 
10 -5 — we will see that this approximation is in fact 
fairly accurate. 

For our numerical calculations, we evolved a scale 
invariant Harrison- Zeldovich spectrum up to redshift 
z = 100 (before which the GCG is indistinguishable 
from ACDM) with CMBfast [19] to correctly include all 
the relevant effects (early super-horizon evolution, pre- 
recombination acoustic oscillations, Silk-damping, etc.), 
with cosmological parameters given by the concordance 
model of [2]. We then used equation (13) to evolve the 
fluctuations from z — 100 until today. Results for a sam- 
ple of a- values are plotted in Figure 1, and show how 
tiny non-zero values of a result in large changes on small 
scales as expected. 

We constrain a by making a x 2 fit of the theoretically 
predicted power spectrum against that observed with the 
2dF 100k Galaxy Redshift Survey [20] as analyzed by 
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Even before solving this, it s obvious that a non-zero 
sound speed, if present for a sufficiently long time-span, 
is going to have a dramatic effect on the fc-dependence 
of the perturbation growth. If c 2 > 0, then fluctuations 
with wavelength below the Jeans scale Aj = ^ n\c 2 \ / G p j= 
will be pressure-supported and oscillate rather than grow. 
This oscillation is confirmed by the numerical solu- 
tions, and is analogous to the acoustic oscillations in 
the photon-baryon fluid in the pre-decoupling epoch. If 
c 2 < 0, corresponding to negative a, fluctuations below 
this wavelength will be violently unstable and grow ex- 
ponentially [18]. 

A key point which has apparently been overlooked in 
prior work is that whereas all the other terms in equa- 
tion (13) are of order unity or smaller, the sound speed 
term (kc s /aH) 2 can be much larger even if the sound 
speed is tiny |c s | <C 1. This is because c s is multiplied 
by the prefactor k/aH which can be enormous, since it 
it the Horizon scale divided by the perturbation scale. 
Defining a critical wavelength A c by 
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FIG. 1. UDM solution for perturbations as function of 
wavenumber, k. From top to bottom, the curves are GCG models 
with a = -10~ 4 , -lCT 5 , (ACDM), 10" 5 and 10" 4 , respectively. 
The data points are the power spectrum of the 2df galaxy redshift 
survey. 



[21]. For each a, we use the best fitting power normal- 
ization to ensure that our constraints come only from 
the shape of the power spectrum, not from the overall 
amplitude which involves mass-to- light bias. To be con- 
servative and stay firmly in the linear regime, we discard 
data with k > 0.3/i/Mpc. We run our code for a fine 
grid of models with — 1 < a < 1 to find the correspond- 
ing x 2 values. The likelihood function e _Ax I 2 is plotted 
in Figure 2. It predictably peaks around a w 0, and 
the observed skewness is simply due to the fact that the 
oscillating solution (a > 0)is easier to fit than the expo- 
nentially unstable solution (a < 0) . Setting Ax 2 = 1 cut- 
off as in a crude Bayesian analysis gives the constraints 
-0.00000081 < a < 0.0000079. 

To place this result in context, Figure 3 shows the 0th 
order constraints from Makler et al. [13] with our new 
strict constraints superimposed. 



IV. CONCLUSIONS 

Above we showed that GCG models with |a| 3> 10~ 5 
are ruled out by observation, since they cause fluctua- 
tions or blowup in the matter power spectrum that are 
not observed. Let us now examine the assumptions that 
went into this and the broader implications. 

First of all, our extremely tight constraints imply that 
that the narrow range of allowed GCG models will be 
completely indistinguishable from ACDM to both to 
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FIG. 2. The likelihood function e~ AX//2 as a function of 
the GCG parameter a. It is sharply peaked around a = 
which is equivalent to the ACDM model. From top to bottom, 
the horizontal dashed lines correspond to A% 2 = 1 and 4, 
respectively. 
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FIG. 3. The graph is showing constraints from previous 
work by Makler el at. Our new constraints from first order 
perturbation theory are superimposed on the plot as shown 



0th order and at the early times when primary CMB 
anisotropics are produced. This means that the corre- 
sponding standard constraints on cosmological parame- 
ters from CMB, SN la etc. apply also to the GCG models 
making the identification fl* n = Q m , so that there are no 
interesting degeneracies between a and other parameters 
that can significantly widen the allowed cv-range. We 
have therefore used standard constraints 0.2 < Q m < 0.4 
for the allowed region in Figure 3. 

Second, our limiting our constraints to linear scales 
k < 0.3ft./Mpc was probably overly conservative. As re- 
viewed in [22,23], there are quite strong constraints on 
the linear power spectrum on much smaller scales from 
weak lensing, from the Lyman a forest and perhaps even 
from lensing of halo substructure [24] which if used would 
tighten our upper limit on a to 10~ 6 , 10~ 7 and 10~ 10 , 
respectively. 

Third, we saw that all that really mattered in equa- 
tion (13) as far as the constraints were concerned was the 
pressure term (A c fc) 2 . This means that our results apply 
more generally than merely to the GCG case: any unified 
dark matter model where p is a unique function of p is 
ruled out if the sound speed is non-negligible, i.e., if the 
function p(p) departs substantially from a constant over 
the range where pressure has an effect — quantitatively, 
if Idlnp/dln p\ > 10~ 5 , again rendering it indistinguish- 
able from standard ACDM. 

In contrast, standard quintessence models have no 
such problems. Although they typically have high sound 
speeds causing oscillations as above, this does not pre- 
vent the dark matter from clustering since it is a separate 
dynamic component. Quintessence models would fail as 
above if there the two components were tightly coupled, 
and this is effectively what happens with UDM since the 
two are one and the same substance. 

To salvage the UDM idea in some form, its pressure 
must not be uniquely determined by its density — not 
even on subhorizon scales. As worked out in detail by 
Hu [18], the effective sound speed can under some cir- 
cumstances differ from the adiabatic sound speed, and 
it is only the former that must approximately vanish to 
satisfy our constraints. Although it may be possible to 
concoct such models, say by introducing another physi- 
cal field upon which p depends and making it fluctuate 
in such a way as to cancel the problematic pressure gra- 
dients, this would be giving up much of the elegance and 
simplicity that gave the unified dark matter idea its ap- 
peal, essentially substituting one extra field for another. 

In conclusion, precision data is gradually allowing us to 
test rather than assume the physics underlying modern 
cosmology. We have taken a step in this direction by 
ruling out a broad class of so-called unified dark matter 
models. Our results indicate that dark energy is either 
indistinguishable from a pure cosmological constant or a 
separate component from the dark matter with a life of 
its own. 
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